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ABSTRACT 


A study was conducted of modeling changes to the XR-3 
Captured Air Bubble Surface Effect Ship digital computer Loads 
and Motion program. These changes included the addition of 
pressure rate and volume rate to the existing 6 degrees of 
freedom equations. Additional equations were developed to 
Simulate the application of a nonpermeable membrane to the 
Plenum of the XR-3 test craft. The objectives of this study 
were to determine whether the addition of pressure rate and 
volume rate equations would improve computer execution time 
and to test some simplified models of the plenum membrane. 
Computer timing improvements were demonstrated and membrane 


modeling results are presented. 
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I. INTRODUCTION 


A. BACKGROUND 

The Captured Air Bubble Surface Effect Ship (CAB SES) 
offers a dynamic new approach to improve surface ship per- 
formance. As with any new approach, extensive testing at the 
design stage of development is required. The U.S. Navy has 
been conducting sea trials on two 100 ton Captured Air Bubble 
(CAB) test craft, the 100-A and the 100-B / Ref. 1.7. In 
addition to scale model testing, another approach to design 
evaluation involves the use of high speed digital computers 
to simulate the dynamic performance of the craft. 

A digital computer Loads and Motion (L&M) program for the 
CAB SES was developed under government contract by the 
Oceanics Corp. / Ref. 2_/7. The L&M program installed at the 
Naval Postgraduate School, W. R. Church, computer facility 
on an IBM 360/67 computer was for a 100 ton displacement test 
craft, the 100-B. Leo and Boncal / Ref. 3_/ modified this 
computer program to represent a smaller 3 ton craft, the XR-3. 
The XR-3 is an operational test craft currently maintained 
and operated by the Naval Postgraduate School / Ref. 4 /. 

Various modifications and changes to Leo and Boncal's 
basic work have been implemented. Finley, Forbes, and Menzel, 
in Refs. 5-7, provides program changes to obtain better 
representation of the Fan maps, bow and stern seal dynamics 


and pitch and roll damping. 
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One adverse characteristic of both the 100-B computer 
program and the XR-3 computer program was the extensive use 
of the digital computer time when operating in simulated sea 
states. Since the 100-B and the XR-3 programs are similiar, 
the computer time analysis provided by Mitchell / Ref. 8 7 
for the 100-B also applies to the XR-3 program. 

Another aid to design is the testing done in towing tanks 
in which scale models are used. Towing tests were conducted 
at NSRDC using a rubber membrane installed on the forward 
portion of a CAB SES model in order to better scale the 
pressure-volume relationship / Ref. 9_/. The addition of 
this membrane improved the pitch and heave characteristics 


of the scale model in sea state operation. 
B. OBJECTIVES 


It is the objective of this work to examine the L&M 
digital computer program and make any changes that would 
improve the program with particular attention to changes 
that would reduce the amount of time required to execute 
the program. 

A second objective is to model the plenum membrane, add 
the required program changes and study the effects of the 


membrane on simulated craft performance. 
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II. PROGRAMMING CHANGES 
A. RATE EQUATIONS FOR PRESSURE AND VOLUME 


ne Sn troduction 

An examination of the dynamics of the XR-3 L&M pro- 
gram indicate that the key parameters of this type of ship 
design are those which describe the dynamic behavior of the 
bubble of the air in the plenum. This pressure of air provides 
approximately 75% of the lift force when the craft is "on the 
bubble." Small changes in bubble volume and pressure there- 
fore cause significant changes in total ship response. These 
variables were examined in an effort to decrease the digital 
computer simulation execution time. 

2. SUBROUTINES RHS and INTGRL 

The force and moment equations for the craft are the 
six degrees of freedom equations plus four auxillary equations 
contained in SUBROUTINE RHS / Ref. 2_/. 

The integration of these equations is performed in 
SUBROUTINE INTGRL. INTGRL uses the Runga-Knutta-Merson 
numerical integration technique with an automatic variable 
step size. If the error tolerances specified are exceeded 
the step size is reduced and the step repeated until all 
error tolerances are met or the step size is reduced to 
i x 107° sec., at which point the program stops. Thus the 
integration values are forced to converge to their correct 


values. All ten equations are integrated serially at the 
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same time step. When all integration error tolerances are 
met INTGRL returns the integrated values and increases the 
Simulation time one time-step increment. Since pressure and 
volume are key variables, integrating their rates should cause 
the system of equations to converge more rapidly to their true 
values and avoid any problems associated with the feedback 
loops in the pressure-volume calculation / Ref. 10_/. 

3. Formulation of Volume Rate and Pressure Rate 

The bubble pressure equation is for an adiabatic 


process / Ref. 2//. 


P =P b (1) 


where the subscripts a and b represent atmosphere and plenum 


values respectively 


P = pressure absolute 
p = density 1. (2) 
y = adiabatic exponent for air 
Supstituting (2) into (4) 
Y 


b a Vo 
a 


Differentiating (3) with respect to time where P, and p, 


are assumed to be constant yields 


. P = . ° 
P = a M Y VM - MV 
ge a ileas = (4) 
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The volume equation for the bubble is computed in RHS in four 


parts and summed together 


V = Vnom - D - WA + WT (7) 


The four terms that determine volume are 


(1) Vnom = the empty plenum volume 


G2), oD = draft term-change in volume due to draft 

(3) WA = Wave term-change in volume due to the presence 
of waves 

(4) WT $= WATSLP term-- correction term for the added 


volume wedge which is a function of speed 


L fet. 6. 


Differentiating equation (4) 


V = -D - WA (8) 
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where 


WA = rate of change of volume due to waves 

D = rate of change of volume due to displacement 
Wr = 0.0 

Vnom =/050 


WA is computed in SUBROUTINE WAVES by taking the difference 

in change of volume due to waves at successive time steps 

and dividing by the time step. Since WATSLP is a function 

of speed, the differentiated value of WATSLP, WT, isa 

function of acceleration. Forward accelerations are normally 
small; therefore, the differentiated value of WATSLP was 
approximated to be zero. Equations (4) and (5) were programmed 


as follows: 


VALUE (11) = -(((XL*WIDTH) -ABW)*. 5*W)-DVWDOT 
VALUE (12) = GAM*VAL(13)*(VALUE(10) /VAL(11)-VALUE(11)/ 
VAL(12) ) 
where 
DVWDOT = rate of change of volume due to waves 


VALUE(10) = Bubble mass flow rate 
VALUE(11) = Volume rate 


VALUE(12) = Pressure rate 


VAL(11) = Bubble Mass 

VAL(12) = Bubble Volume 

VAL(13) = Bubble Pressure 

GAM = Adiabatic Exponent for air 
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B. ELIMINATION OF SUBROUTINE DMINV 

In an effort to improve program timing, the suggestions 
given in Ref. 7 were undertaken. One recommendation suggested 
the removal of SUBROUTINE DMINV in the 100B L&M program by 
making changes to SUBROUTINES INCON and RHS. An examination 
of the L&M program indicated that these changes also could 
be applied to the XR-3 program. Consequently the following 


was deleted: 


2a2 veDOe Zi IIIS! 156 
DO 221 N = 1,6 
21 6=— AC LGN) =" 0,0 
DOs 213.N = 19 


213 , A(N,N) ,= AM 
ACE W) = ALKE 
AG, Siva as nlee 
A(6,6) = AIZZ 
A(4,6) = -AIX2 
A(6,4) = -AIX2 
AIMAX = AMAXI(AM, AIXX,AIYY,AIZZ,ABS(AIXZ) ) 


DO 214 1 = 1,6 
DO 214 J = 1,6 
214 A(I,J) = A(I,J)/AIMAX 
CALL DMINV(A,G,D,) 
DO) 2051 9/=7156 


prsie AGI,d) A(I,J) /AIMAX 
IF) (DsNEJO.0") Go" TO To 
WRITE (6,216) 

STOP 
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In place of the above the following was added: 


215  AMASSI = 1.0/AM 
D = 1.0/( AIXK* AIZ2-AIXZ* AIXZ) 
DIXX = AIXX*D 
DIXZ = AI2ZZ*D 
DIZZ = AIZZ*D 
AIYYI = 1.0/ATYY 
GO TO 10 


Linkage between INCON and RHS was provided by the following 


COMMON /ATRIX/AMASSI , AITYYI, DIXX, DIXZ, DI2Z 


In SUBROUTINE RHS the six element matrix GF(J) was deleted 
and the following identifiers were substituted for the sum- 
mation of forces: SUMX, SUMY, SUMZ, SUMK, SUMM, SUMN. In 


addition, the following deletion was made. 


DOr 1° r= 1,6 

VALUE(I) = 0.0 

DOr d=" 1,6 

VALUE (I) = VALUE (I) + A(I,J)*GF(J) 
1 CONTINUE 


Substituted for the above DO LOOPS was the following 


VALUE (1) = SUMX*AMASSI 

VALUE (2) = SUMY*AMASSI -R*U 
VALUE (3) = SUMZ*AMASSI 

VALUE (4) = SUMK*DIZZ + SUMN*DIXZ 


Ib 


VALUE (5) 


SUMM* ATYYI 
VALUE (6) 


SUMN* DIXX + SUMK* DIXZ 
These changes eliminated SUBROUTINE DMINV and several DO LOOPS. 


C. EXPANSION OF WAVE COMPONENTS 


Since the program was being tested with sea state, the 
method used for the wave generation was studied. SUBROUTINE 
INCON provides a means of introducing sea-state by individual 
wave components and amplitudes. In addition, the subroutine 
will accept an average height and lowest and highest wave 
frequency and/or wave period and generate up to 10 wave 
components. 

These wave components represent a sampling of the spectural 
energy density of the given sea state condition. Reference 
11 suggests’ for irregular sea, a minimum of 15 - 20 com- 
ponents are required to simulate irregular sea conditions. 
The program dimension> statements were accordingly increased 
from 10-20 to determine if any significant changes could be 


observed by increasing the number of wave components. 


D. MISCELLANEOUS 


Each subroutine was examined with regard to efficiency 
of coding and changes were implemented such as multiplication 
in lieu of the use of the exponential. For example, IBM 360/67 
fixed point exponentiation requires that the natural logarithm 


be generated and then a logarithm be computed, as compared 
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to straight forward multiplication. When such exponentiation 
is nested in DO LOOPS, a significant time savings can be 


realized by the use of multiplication. 


We 


III. TEST PRODECURES 


A. COMPUTER SYSTEM 


The digital computer used throughout this study is an 
IBM 360/67 model, VERSION I, located at the NPS W. R. Church 
Computer Facility. The system's hardware configuration pro- 
vides 2M bytes of core. This main core is composed of IBM 
core storage devices (IBM 2365 Mod. 12) and compatable but 
slightly different Lockheed core storage devices (MM 365). 
The Lockheed core devices are about 18% slower in execution 
time than IBM core devices. Since the system is time-shared 
and the core is contiguous, there is no practical way to 
determine in which portion of the main core the program is 
residing. Therefore, exact timing of a program is not easily 
achieved. In this report all timing values must be interpreted 
with this anomaly in mind. Among the options available under 
this operating system are the FORTRAM G and FORTRAN H compilers 
7 Ref. 12_/7. Mitchell reported substantial time savings using 
the H compiler / Ref. 8_/7. Since one of the main objectives 
was to reduce execution time, computer model testing was done 
using both compiler options. Runs made under different 
compilers will be noted. 

A word of caution on the use of the H compiler. The H 
compiler produces its timing improvements and core reduction 


by coding optimization. This means the compiler re-arranges 
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the Fortran code. This re-arrangement may move computations 
outside of DO LOOPS for example. Applicable portions of 


Ref. 13 should be read and understood before using this feature. 


B. TEST CONDITIONS 


The L&M program for the XR-3 used in this report is the 
one given in Ref. 7. This program was then modified to 
incorporate the changes previously described. The basic 
program / Ref. 7_/ computes pressure using eq. (1) and volume 
using eq. (7) and is referred to as. PROGRAM ONE. PROGRAM ONE 
was used as reference for comparative purposes. The addition 
of the pressure rate and volume rate equations (6) and (8) and 
the elimination of SUBROUTINE DMINV constituted major modifi- 
cations to the basic program. This version of the Loads and 
Motion program was called PROGRAM TWO. PROGRAM THREE was 
PROGRAM TWO compiled and executed with the FORTRAN H compiler. 

Five test conditions for a given speed were established 
to provide a full range of disturbances for testing the L&M 
programs. These test conditions are tabulated in Table I. 

In Condition One small disturbances were introduced by 
changing the initial conditions of draft (DS) and pitch (0). 
These changes consisted of decreasing the draft and pitch 
from previously computed steady state values by the magnitude 
shown in Table I. In addition, the rudder was held at zero 
and the sea was calm. In Condition Two, the perturbations 


were the same as Condition One except that now the rudder 


al 
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TABLE I 


TEST PERTURBATIONS 


ADS AQ RUDDER SEA STATE 
Condition One GO. 1 “inen rOg 0.0 Calm water 
Condition Two 0.1 03 20° Calm water 
Condition Three 0.0 0.0 0.0 Regular One 
Condition Four Oe) 0.0 0.0 Regular Two 
Condition Five 0.0 0.0 0.0 Irregular One 


was displaced to the right 20 degrees. The rudder was offset 
at the start of the simulation and held through the run as 
this value. In Condition Three, Four and Five the only dis- 
turbances used were the sea state conditions. 

Sea state is generated in SUBROUTINE WAVES. SUBROUTINE 
INCON generates wave components and amplitudes based on 
several different type input parameters. ZL See users manual 
her, 7./. Regular sea one and two are simple sinusoidal 
waves. Irregular sea state conditions are simulated by the 
addition of several regular sinusoidal wave components. Sea 
state used throughout this study is true sea state and is 
not scaled in any manner. 

The L&M program provides two options for propulsion, 


either constant speed or constant thrust. The constant 


ee 


thrust option was selected because this is the normal oper- 
ating mode for the XR-3 test craft. 

Integrator tolerance levels were maintained at the values 
given in Ref. 7. However two new error tolerances had to be 
determined for the integration of the pressure and volume 
rates. These new levels were determined by first selecting 
a very tight value .000001 which was increased in increments 
until a change in output could be noticed. The value was then 
decreased by a factor of 10. Pressure and volume tolerance 
levels were then set at this value which was found to be .0001. 

All other variables such as draft.and pressure were ini- 
tialized from the steady state undisturbed conditions for that 
specific starting speed. The steady state conditions were 
evaluated by first using the constant speed propulsion option 
and the initial values in Ref. 7. The program was then run 
until the key input variables - pitch, draft, thrust and 
pressure - reached steady state. These values were then 
read in and propulsion was switched to constant thrust. 

The program was run again to insure that steady state con- 
ditions existed. These steady state values are shown in 


Table II. 
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TABLE II 


STEADY STATE CONDITIONS 


Speed(Knots) 10 20 30 
Pressure 23.93 24,84 24.84 
Dratt 857 6.12 5.34 
Thrust (each engine) 200.31 21:8. 1.7 287.22 
Eaten A6z 48 . 26 
Pressure in psfg Thrust an £4. Lb. 

Draft in inches Pitch in degrees 
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IV. MEMBRANE ADDITION 


A. INTRODUCTION 


During towing tank tests conducted at NSRDC, a rubber 
membrane was installed on the forward bow of the XR-1 scale 
model / Ref. 9_/. The installation of this membrane improved 
the pitch and heave characteristics of the scale model during 
simulated sea conditions. Reference 9 provided the following 
emperical expression which relates the Pressure and Volume 


of the membrane. 
Vr=(KP (9) 


K was experimentally determined and has the value .317 when 


scaled to the XR-3 test craft dimensions. (See Appendix A). 


B. MATHEMATICAL MODELING 


The approach used to include the effects of this membrane 
was to develop mathematical expressions which could be added 
to the already existing pressure, volume and air mass equations 
in the L&M program. Five mathematical expressions were devel- 
oped, tested and referred to as Models One, Two, Three, Four 
and Five. They are explained below and the results discussed 
in Section V. 

1. Model One 

Model One is based on the premise that by correcting 


the plenum volume with the membrane volume, the equations 
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in SUBROUTINE RHS would show the effect of the membrane. 
Model One added the following equation to PROGRAM ONE; 


Ve =e 2 (10) 


Then using (10) to correct the volume of the plenum and 


assuming a = Ps 


V=V.+V ‘@alb) 


where V, is given by eq. (7). 


2. Model Two 
Model Two is an extension of Model One. Based on the 
premise that by correcting the rate of change of plenum volume 
with the rate of change of membrane volume, the rate equations 
given in SUBROUTINE RHS and subsequent integration of these 
terms would show the effect of the membrane. Model Two added 


the following rate equation to PROGRAM TWO. 
{= 2 (12) 


Then using (12) to correct the rate of change of plenum volume 
and assuming 
Pa . Fo 


V =Vi+V 


where Vy is given by eq. (8). 
3. Model Three 
Model Three was based on the premise that the volume 
rate eq. (8) and the pressure rate eq. (6) as computed without 
the membrane could be corrected for the membrane by adjusting 


these terms by appropriate factors within PROGRAM TWO, 
26 


Correction terms were sought that would be a function of 
the membrane. A logical candidate for this correction factor 
for the pressure was the ratio of rates of change of the 
membrane volume to plenum volume multiplied by plenum pressure 
rate. 

V 
Por == PL (13) 
V 
he) 
where PCT is the correction factor to be added to the plenum 
pressure rate term P. before integration. It was assumed 


that the pressure and pressure rate of the plenum and membrane 


were equal. 


P= P, (14) 
Pee 5) 


It was also assumed that the membrane volume rate and pressure 


rate were related by 


Me = KP (16) 


Substituting (16) into (13) and using (15) 
Pa 
PCT = K —~—_ (17) 
ue 
The volume rate correction term was assumed to be the 
ratio of the rate of change of membrane pressure to the rate 


of change of plenum pressure times the rate of change of mem- 


brane volume. 
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VCT = ete Vn (18) 
b 


where VCT is the correction factor to be added to the plenum 


volume rate term V5 before integration 


Since ee = P. and Ve = KP. 


Then from (18) 


SON eee (19) 


Then using equations (17) and (19) the rate equations are 


V 


v, + VCT (20) 


P 


P+ PCT (21) 
where V, is given by eq. (8) and P, by eq. (6). 


Consider the following. Since K is a function of the 
elasticity of the membrane, if K = 0.0. Equations (20) and 
(21) reduce to the no membrane case. 

If the membrane expands at the same rate that the 
volume of the plenum contracts the system volume rates is 
zero and the system pressure rate is also zero. 

4, Model Four 

Again using PROGRAM TWO, Model Four was based on the 
premise that the mass of the membrane air and mass flow rate 
of membrane air would be important factors in any correction 


term to be applied to plenum volume rate and pressure rate. 
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Again it was assumed that the membrane effect could 
be described as an adiabatic reversible process. Then using 


Equation (6) applied to the membrane, 


m 


M V 
m m 


M v 
m -_ 


ca ( 22) 


where m denotes membrane 


M._ = Mass flow rate of the membrane air 


M = Mass of the membrane air 


Again it was assumed that the pressure and pressure rate of 


the membrane and plenum were equal. 
P =P (23) 
P =P ( 24) 


It was also assumed that the membrane volume rate and pressure 


rate were related by 


aes Be (25) 


Equation (25) was obtained by differentiating 
V =KP (26) 
Substituting (25) and (26) into (22) yields 


M ee P ae 
M ) P: 
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One term is not available, Ma Two possibilities 
were considered to determine this value. The first possibil- 
ity would be to integrate M and use the result of the inte- 
gration from the preceding time step. The other possibility 
would be to hold temperatures constant in the membrane, and 
use the ideal gas law to determine Mn for varying pressure 
and volume. This last condition would be in error; however, 
as a first approximation, it would produce acceptable results. 

Because of the past history of difficulty with inte- 
grator 10, plenum air mass rate, the second possibility was 
chosen. Evaluation of the terms in the ideal gas law for 
a temperature of 68° degrees excluding mass, pressure and 


volume produced the constant. 
4.g X 107’ slugs/psf - ft? 


therefore; 


M = (4.9% 1077) (P_)(v)) ( 28) 


Equations (22), (25) and (27) describe the membrane. 
These equations are used to correct the mass rate and volume 


rate of the plenum without the membrane. 


: 7 ae 
M. M+ M, (29) 
yo ee we. (30) 

Ss mh 

where 
s = denotes the system of plenum plus membrane 
M = mass flow rate of plenum 
V = volume rate of change of plenum 
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Using (29) and (30) in equation (6) the corrected plenum 
pressure rate is 

; M V 

Te Goi sine ea 
where PL. M and V are values obtained from the previous time- 
step calculation. 

5. Model Five 

Model Five is the application of equation (10) to 

PROGRAM TWO. It is based on the premise that the volume of 


the membrane could be introduced into the plenum pressure rate 


equation, and represent the membrane effect. 


Vv 


HT} 
A 
td 


(31) 


m m 


il 
< 
+ 
< 


u b m (32) 


Equation (31) and (32) are computed as in Model One. Then 


equation (32) is used in the pressure rate equation. 


: M V 
Ewe 0 lala Ca (33) 


Be 


lL. > 
7 (2) pet¢supe ad : (GE? 


V. DISCUSSION OF RESULTS 
A. PRESSURE RATE AND VOLUME RATE CHANGE 


The addition of the pressure rate and volume rate equa- 
tions was implemented without any major difficulty. In order 
to provide a basis of comparison for computer time analysis 
a quality factor "Q" was established. "Q" is defined to be 
the ratio of CPU execution time to problem time. For example, 
if the problem specified 1 minute of simulation and the CPU 
execution time (GO STEP in IBM COMPUTER SYSTEMS) was 5 minutes 
Q = 5.0. PROGRAMS ONE, TWO and THREE were tested with identi- 
cal input parameters. The results are summarized in Table III 
and Table IV. 

Runs were compared using PROGRAM ONE as a reference. No 
differences in computed output were noted for Condition One 
thru Four and only small changes in computed values noted with 
Condition Five. The integrator error tolerance levels were 
maintained constant for all runs. 

As can be seen from the tabulated values, the addition of 
pressure rate and volume rate equations reduces the CPU execu- 
tion time considerably for small disturbances. However, as 
the sea state is increased the CPU time increases. Apparently 
other factors in the program begin to dominate the total execu- 


won time. 
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TABLE III 


TIMING QUALITY FACTORS 


20 Knots 
Program One Program Two Program Three 

Condition One 25 4,8 4,50 
Condition Two 13'.35 6n6 5.50 
Condition Three Page 8) 2On Lye 
Condition Four 42.8 28.01 2250 
Condition Five TAA en 9) Sismle 54.0 

TABLE IV 


TIMING QUALITY FACTORS 


30 Knots 
Program One Program Two Program Three 
Condition One 190.6 Te5 4.5 
Condition Two 18.0 14.25 Si 
Condition Three 24.0 22.30 21.0 
Condition Four 45.0 30.0 23..0 
Condition Five 75.0 59.2 be 0 


25) 


This supports the results of Reference (7) that for sea 
conditions the bulk of computation time is apparently in 
the SUBROUTINE WAVES. 

It was observed that at irregular sea state two and 
regular sea state three water came in contact with the top 
of the plenum. The L&M program does not include the effect 
of this phenomenon in the program simulation; therefore, 
testing was not conducted at higher sea state conditions. 
The fact that water hits the top of the plenum is a reasonable 
condition under certain circumstances. For example, if the 
pitch angle is 3 degrees up the stern sinks over 6 inches. If 
tne draft of the craft at that time is 11 inches and if the 
preceding wave passing the stern station is 5 inches, water 
should hit the top of the plenum at the stern station. The 
effect of this condition on total craft behavior in terms of 
forces and moments is not known. The effect of water contact 
would depend upon such factors as the length of time the 
wave is in contact with the top of the plenum and the area 


of the plenum in contact with the wave. 
B. INCREASING WAVE COMPONENTS 


Tests were made for irregular seas one and two using 8, 
10, 15, and 20 wave components. No apparent change in output 
could be determined. The time at which changes in output 
occurred changed, but no significant reduction in output 
values were noted; that is the time when wave components were 


in-phase or out-of-phase was different, but the craft response 


34 


was the same. It was felt that if more components were used, 
the energy (in this case, the magnitude and frequency of the 
disturbing force) would be distributed over many components 
and that the irregular combination of these components would 
spread the energy of the disturbance over a larger time base. 
This spreading of the energy would then reduce the magnitude 
of the disturbance and that in turn would provide less rapid 
rates of change. With less rapid changes in the disturbances, 
the integration would take less time. The over-all time increase 
due to more wave components would be offset by the reduction 
in integration time. This is not the case; using more wave 
components increased CPU time. Whether 10 or less wave com- 
ponents are sufficient to accurately represent irregular seas 


was not tested. 
C. MEMBRANE MODELS 


A large amount of time was expended in attempting to provide 
adquate modeling of the plenum membrane. Model One using 
Fe =K P. in PROGRAM ONE did not work at all. The program 
would begin to execute and then stop when the integration step 
size was less than 1 X 107°, 

As an aid to find the reason for this failure a computa- 
tional flow diagram was drawn (See Figure 1). The addition 
of the correction term Ves creates an algebraic inner loop as 
shown by the dashed line in Fig. 1. Using a small perturbation 


of D in calm water and using the print switches in RHS as well 


as using the DEBUG option available with the G compiler, it 
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was observed that with Model One the values of pressure were 


divergent. It is not known why pressure values diverge; 


however, it was noted that at certain times the values of 
pressure and volume were increasing or decreasing together, 
that is when pressure was increasing while volume was increas- 
ing. 

The computational flow diagram for Model Two (Figure 2) 
Shows the inner connections of the pressure, volume and mass 
variables. Model Two exhibited the same symptoms and the same 
results as Model One, except in some cases pressure rate and 
volume rate rose or fell together at the same time step. 

The failures of Model One and Two led to Model Three. It 
was felt that by correcting both pressure rate and volume rate 
the divergence exhibited in Models One and Two could be correc- 
ted. Model Three also failed; however, it provided the motiva- 
tion to consider mass rate into the membrane. The computations 
performed in Model Three is shown in Figure 3. 

Model Four was marginally successful. The program exe- 
cutes and slightly reduces the center of gravity acceleration. 
the reduction of the center of gravity acceleration was chosen 
as the criterion of whether the model was representing the mem- 
brane. This model was tested in sea state one and sea state 
two and provides slight reduction in the maximum values in 
the center of gravity acceleration. It did not reduce all of 


the center of gravity values. The computational flow diagram 
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of Model Four is given in Figure 4 and shows no algebraic 
loops in the RHS computations of pressure and volume rates. 

Model Five was the most successful model. The computa- 
tional flow is shown in Figure 5. Model Five reduced the 
center of gravity acceleration and for small disturbances, 
the pitch angle and pitch rate. For sea state one a signifi- 
cant reduction was observed. (See Figures 6 and 7). However, 
for sea state two the center of gravity accelerations were 
greater for the membrane when compared to the no membrane 
case. (See Figures 8 and 9). No reason for this could be 
found. The only differences in these two runs was the sea 
state. 

One possible reason for the sea state two results, though 
not yet proven, could be a resonance effect at the sea state 
encounter frequency. Concurrent work in this area / Ref. 14/7 
investigates the conditions for the existence of any resonance 
phenomena. Preliminary investigations indicate that ata 
wave encounter frequency of approximately 4 rad / second 
resonance in pitch angle occurs. The encounter frequency for 
regular sea state two is 5 rad / second. Whether resonance 


phenomena could cause this effect, however, was not determined. 
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FIGURE 8 
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FIGURES 9 


X-SCALE=5.00E+00 UNITS INCH. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. PROGRAM TIMING 


The inclusion of pressure rate and volume rate equations 
improved the execution time of the L&M program without decreas- 
ing the program accuracy. The changes to the program were not 
difficult to implement. Further effort to improve execution 
time of the program should be directed toward changes to 
SUBROUTINE WAVES. In connection with this, a study of the 
number of required sidewall and bow and stern seal stations 
should also be conducted. By reducing the number of sidewall 
stations the execution time should decrease. A study of the 
losses in accuracy, however, would be necessary. 

Another facet for investigation within SUBROUTINE WAVES 
would be a study to determine if the trignometric functions 
could be represented by the in-line programming of the expan- 
sion of sine and cosine functions, vice the use of computer 
library functions. Since some sine and cosine functions are 
nested in DO LOOPS small savings in time would accumulate. 
Again, the losses in accuracy to this change would have to be 
determined. 

The elimination of SUBROUTINE DMINV did not produce any 
noticeable time savings but did save approximately 2 K bytes 


of core. 


4? 


sy@ irae 
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The expansion of the number of wave components did not 
demonstrate any time savings; however the optimum number of 
wave components required to simulate irregular sea conditions 
should be determined. A possible course of action would be 
to generate wave components, calculate the magnitude and 
frequency and compute the energy spectral density. This 
energy density spectrum could then be compared to the Pierson 
and Neumann or Moskowitz model [ Ref. smear /A By determin- 
ing the optimum number of wave components, a savings in CPU 


execution time could result. 
B. ADDITIONAL MODELING 


The fact that waves can hit the top of the plenum was 
discussed in Section V. Whether this effect produces signif- 
icant forces and moments is unknown. Since the time that the 
wave is in contact with the plenum and the location (station 
numbers) is known and since the area of the plenum in contact 
with the waves could be approximated, (no greater than the 
difference between wetted stations) additional correction 
terms could be developed to determine whether significant 
forces and moments are generated by wave contact. The drag 


of the wetted plenum top could also be considered. 


C. MEMBRANE MODELING 


Of the five membrane models developed in this study, 
Model Five is the most usable. The accuracy of this model 


must be determined and the model must be validated. Whether 
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mass and mass flow rate terms could improve this model also 
needs to be determined. In this regard a detailed study of 
adiabatic-reversible unsteady flow conditions would be useful. 
This would mean the establishment of a suitable control volume 
and the solution of energy-work equations. 

The membrane modeling in this study was for a continuous 
non-permeable membrane. Additional modeling for the case of 
a hole cut into the membrane should be attempted. The equation 
for a sharp-edge orifice could be used to determine air loss. 
The area (size) of the hole can be determined by considering 
as a first approximation the ratio of surface area of a hemi- 
sphere to volume of the hemisphere. As an alternative, the 
experimental results of the towing tank tests conducted at 
NSRDC, on the XR-1, using suitable scale factors, could be 
used to determine empirical correction factors to be applied 
to the Loads and Motion program. 

The occurrences of algebraic loops and summation of feed 
forward terms as shown in the computational flow diagrams and 
the failure of the models in which they occur indicate that 
further modeling of the membrane is required. To prevent model 
failures, elimination of the effects of the computational loops 


appears to be necessary. 


no 


APPENDIX A 
MEMBRANE SCALE FACTOR 


The empirical relationship between pressure and volume 


for the membrane given in Ref. 9 is for the XR-1 test craft. 


Le = K Pa (Al) 
where 
_ 40 in? 
a Psf 


It was therefore necessary to scale the constant K for the 
XR-3 test craft. 

In the following development the subscript 1 refers to 
the XR-1 and subscript 3 refers to XR-3. 

Since XR-1 and XR-3 are similar it was assumed that all 


dimensions scale by the ratio of plenum lengths 


oy (A2) 
“ 
where 
L, = 3 
3 20 £t 
L = 5,4 2b, 


It was assumed that the ratio of mass to volume was constant 


M M 
sed Eee ae 
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mae ratio of craft volume is 


where g is the gravitational constant. 


Substituting (A2) and (A3) into (45) 


P v4 
a7 ee e = =m 
1 te 


Using (A4) and (A6) in (Al) 


3 
ee ib a eee 
ab 1 Po/ Na 
“i 2 
- = } K 


Converting in.’ to ft.’ and using (A2) 


2 
K 7 a) ee tt?/psf 


ppl 


(A4) 


(A5) 


(A6) 


(A7) 


(A8 ) 


(A9 ) 


APPENDIX B 


LISTING OF FORTRAN XR-3 


L&M COMPUTER PROGRAM 


PROGRAM TWO MODEL FIVS 
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